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ABSTRACT 

An infinite family of axially symmetric thin disks of finite radius is presented. The 
family of disks is obtained by means of a method developed by Hunter and contains, 
as its first member, the Kalnajs disk. The surface densities of the disks present a 
maximum at the center of the disk and then decrease smoothly to zero at the edge, 
in such a way that the mass distribution of the higher members of the family is more 
concentrated at the center. The first member of the family have a circular velocity 
proportional to the radius, representing thus a uniformly rotating disk. On the other 
hand, the circular velocities of the other members of the family increases from a value 
of zero at the center of the disks until a maximum and then decreases smoothly until 
a finite value at the edge of the disks, in such a way that for the higher members of 
the family the maximum value of the circular velocity is attained nearest the center 
of the disks. 
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1 INTRODUCTION 

A fact usually assumed in astrophysics is that the main part 
of the mass of a typical spiral galaxy is concentrated in a thin 
disk l|Binnev fc Tremaine (1987. )). Accordingly, the obten- 
tion of the gravitational potential generated by an idealized 
thin disk is a problem of great astrophysical relevance and 
so, through the years, different approaches has been used to 
obtain such kind of thin disk models. 

Wyse and Mayall ()1942i 'l studied thin disks by superpos- 
ing an i nfinite family of elementary disk s of d ifferent radii. 
Brandt fig&o') and Brandt and Belton (jlQed ) constructed 
flat galaxy disks by the flattening of a distribution of mat- 
ter whose surface of equal density were similar spheroids. A 
simple potential-densi ty pa ir for a thin disk model was in- 
trodu c ed by Kuzmin ([1953) and then rederived by Toomre 
(|l963l : ll964l ') as the first member of a generalized family of 
models. 

The Toomre models are obtained by solving the Laplace 
equation in cylindrical coordinates subject to appropriated 
boundary conditions on the disk and at infinity. The Kuzmin 
and Toomre models of thin disks, although they have sur- 
face densities and rotation curves with remarkable proper- 
ties, represent disks of infinite extension and thus they are 
rather poor flat galaxy models. Accordingly, in order to ob- 
tain more realistic models of flat galaxies, is better to con- 
sider methods that permit the obtention of finite thin disk 
models. 
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A simple method to obtain the surface density, the grav- 
itational potential and the rot ation curve of thin disks of 
finite radius was developed by iHunteil (ligeal ). The Hunter 
method is based in the obtention of solutions of Laplace 
equation in terms of oblate spheroidal coordinates, which are 
ideally suited to the study of flat disks of finite extension. By 
superposition of solutions of Laplace equation, expressions 
for the surface density of the disks, the gravitational poten- 
tial and its rotational velocity can be obtained as series of 
elementary functions. 

The simplest example of a thin disk obtained by means 
of the Hunter method is the well known Kalnajs disk 
(|Kalnaisl(|l972h '). which can also be obtained by fla t tening a 
uniformly rotating sphero i d (IWvse fc Mavalll (|l942l ): lBrandtl 
il96(i ): iBrandt fc BeltonI (| 19621 ) ). The Kalnais disk have a 
well behaved surface density and represents a uniformly ro- 
tating disk, so that its circular velocity is proportional to 
the radius, and its stability properties have been extensively 
studied (s ee, for instance. Hunter (fl96 3; 1965), Kalnajg 
(|1972l ) and lKalnais fc Athanassoula-Georgalal (|l97^)). 

In this paper we use the Hunter method in order to ob- 
tain an infinite family of thin disks of finite radius. We par- 
ticularize the Hunter general model by considering a family 
of thin disks with a well behaved surface mass density. We 
will require that the surface density be a monotonically de- 
creasing function of the radius, with a maximum at the cen- 
ter of the disk and vanishing at the edge, in such a way that 
the mass distribution of the higher members of the family 
be more concentrated at the center. 

The paper is organized as follows. In Sec. 2 we present 
a summary of the Hunter method used to obtain the thin 
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disk models of finite radius and also we obtain the general 
expressions for the gravitational potential, the surface den- 
sity and the circular velocity. In the next section, Sec. 3, we 
present the particular family of models obtained by impos- 
ing the required behavior of the surface densities and then, 
in Sec. 4, we analyze its physical behavior. Finally, in Sec. 
5, we summarize our main results. 



2 GENERAL FINITE THIN DISK MODELS 

In order to obtain finite axially symmetric thin disk mod- 
els, we need to find solutions of the Laplace equation that 
represents the outer potential of a thin disklike source. Ac- 
cording with this, we need to solve the Laplace equation for 
an axially symmetric potential. 



H 



0, 



(1) 



where {R, 4>, z) are the usual cylindrical coordinates. We will 
suppose that, besides the axial symmetry, the gravitational 
potential has symmetry of reflection with respect to the 
plane z = 0, 



^{R,z) = ${R,-z) 



(2) 



so that the normal derivative of the potential, d^/dz, sat- 
isfies the relation 



9$ 9$ 



(3) 



in agreement with the attractive character of the gravita- 
tional field. We also assume that d^/dz do not vanishes 
on the plane z — 0, in order to have a thin distribution of 
matter that represents the disk. 

Given a potential z) with the above properties, the 
density E(_R) of the surface d istribution of matter can be 
obtained using the Gauss law (jBinnev fc Tremainj |l983)). 
So, using the equation (|3]), we obtain 



E(i?) = 



2nG 



a7 



(4) 



Now, in order to have a surface density corresponding to a fi- 
nite disklike distribution of matter, we impose the boundary 
conditions 



oz 

oz 



= 0; 



R<a, 



R> a, 



(5a) 
(5b) 



so that the matter distribution is restricted to the disk z = 0, 
0<R<a. 

We introduce now the oblate spheroidal coordinates, 
whose symmetry adapts in a natural way to the geometry of 
the model. This coordinates are related to the usual cyli n- 
drical coordinates by the relation (|Morse fc FesbachI l|l953h ). 



R^a^{l+e){l-V^), 
z = a^-n, 



(6a) 
(6b) 



where < ^ < oo and —1 < rj < 1. The disk has the 
coordinates ^ = 0, < t?^ < 1. On crossing the disk, rj 
changes sign but does not change in absolute value. This 
singular behavior of the coordinate rj implies that an even 



function of ry is a continuous function everywhere but has a 
discontinuous rj derivative at the disk. 

In terms of the oblate spheroidal coordinates, the 
Laplace equation can be written as 



(7) 



and we need to find solutions that be even functions of rj 
and with the boundary conditions 



'i',dO,v)=F{rj), 

■J',.(C,o) = o, 



(8a) 
(8b) 



where F{rj) is an even function which can be expanded in a 
series of Legend re polynomials in the interval — 1 < »7 < 1 
jBatemanl (|l944l ')'). 

According with this, the Newtonian gravitational po- 
tential for the exterior of a finite thin disk wit h an axi- 
ally symmetric matter density can be written as (|Batemanl 
(0944,)), 

oc 

3>(C,r?) =-EC2n?24C)^2„(77), (9) 
11=0 

where C2n are arbitrary constants, P2n{ri) and g2n(C) = 
i'^"+^Q2„[i^) are the usual Legendre polynomials and the 
Legendre functions of second kind, respectively. With this 
general solution for the gravitational potential, the surface 
matter density is given by 



E(i?) 



2iiaGrj 



^C2„(2n+l)g2n + l(0)P2n(?)) (10) 



and, as we will shown later, the arbitrary constants C2n 
must be chosen properly so that the surface density presents 
a physically reasonable behavior. 

Besides the matter density, another quantity commonly 
used to characterize galactic matter distributions is the cir- 
cular velocity V{R), also called the rotation curve, defined as 
the tangential velocity of the stars in circular orbits around 
the center. Now, given "I>(_R, z), we can easily evaluate V 
through the relation 



dR 



in such a way that, by using ((9)1, we obtain 
R^ 



V%R) 



(a2 - i?2)i/; 



^C2ng2n(0)P:^„(77). 



(11) 



(12) 



for the circular velocity. 



3 THE GENERALIZED KALNAJS DISKS 

We will now particularize the above general model by consid- 
ering a family of finite thin disk with a well behaved surface 
mass density. We will require that the surface density will 
be a monotonically decreasing function of the radius, with 
a maximum at the center of the disk and vanishing at the 
edge. In order to do this, we impose the conditions 



E(a) 
E(0) 



0, 



(13) 
(14) 
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and we also require that 

M = 2-K 



T.{R)RdR, 



where M is the total mass of the disk. 

Now, by using the boundary condition 
density can be written in the form 



(15) 



, the surface 



E(7?) 



2-KaGrj ' 



(16) 



where F{ri) is an even function of 77, monotonically increas- 
ing at the interval < r; < 1, and such that 

lim = 0. (17) 

'7->o r] 

Furthermore, we must to impose the condition 

Fir,)dri=^, (18) 

in agreement with (|15p . 

A simple function F(ri) that agrees with a ll the a bove 
requirements was given by iLetelier fc Oliveiral (|l987l ) and 
can be written as 



a 



(19) 



where, in order to fulfill the condition (|17p . we must take 
m > 1. With this particular choice of -F(r?) we obtain an in- 
finite family of finite disks with surface mass densities given 
by 



S,„(i?) 



(2m + 1)M 



1-5 



2-1 m-1/2 



(20) 



As we can easily see, the dis k with m = 1 corresponds to 
the well known Kalnajs disk (|Kalnaisl (|l972l )'). Accordingly, 
this family of finite thin disks can then be considered as a 
generalization of the Kalnajs disk 

Now, from the equation pU|l . the function F{ri) can be 
written as 



with 



F(,7) = ^i^2nP2n(r?), 



K2n = (2n + l)g2n + l(0)C2„. 



(21) 



(22) 



The coefficients K2n are founded, by using the orthogonality 
property of the Legendre polynomials, through the expres- 



K2n 



An + 



F{ri)P2n{v)dv- 



(23) 



The above equation can be expressed as (|Batemanl (Il953t )') 
MG r7r^/2(4n + l)(2m + l)r(2m + 1)" 



K2n = 



2a 



22™r(l + m~ n)r(m + n + 



so that, using the gamma function properties, we obtain: 



C2n = 



MG 
2a 



tn + l)(2m-h 1)! 



22'"(2n + l)(m - ri)!r(m + n+ f )?2n+i(0) 



for n < m and C271 = for n > m. 



4 BEHAVIOR OF THE MODELS 

With the above values of the C2n we can compute the dif- 
ferent physical quantities that characterize the behavior of 
the models. So, for instance, the gravitational potential of 
the first three members of the family are given by 

MG, 



M^,V) = -^[cot-i C + A{3v' - 1)], 

M^,V) = cot 3r?'_l 

a 7 

-I- B(35r?''-30r?^-f3)], 
-l>3(^,r?) = ^[cot i5+-_(3r?2-l) 

+ ^(35r,*-30r,V3) 



(24a) 



(24b) 



-f C(23l77'*- 31577* -f 10577^-5)], (24c) 



where 



^=i[(3e' + l)cot-ie-3C], 



(25a) 



B = ii^I(35e* + 30^' + 3) cot-i i - 35e - f C], (25b) 



231C'' - 238r - -^Cl, 



(25c) 



with similar, but more involved, expressions for greater val- 
ues of m. 

By taking ^ = at the above expressions, we can obtain 
the value of the potential at the plane of the disk. Now, is 
easy to see that, for all the members of the family, we will 
obtain finite values for these quantities. In particular, for 
the first member of the family, the potential at the disk is 
given by 



^ , „ StvMG „2 



(26) 



for R < a, and this expression is completely equivalent to 
the corresponding expression in iKalna^ (|l972l '). 

In the same way, we obtain for the surface mass densities 
of the first three members of the family the expressions 

(1 - -R ) , (27a) 



27ra2 " 

_ ~2s^3/2 
7M n2N5/2 



(27b) 
(27c) 



and for the corresponding circular velocities the expressions 

(28a) 

^2 = fl" (4 - 3/?"), (28b) 

256a - - " (28c) 

where has been introduced the dimensionless radial variable 
R = R/a. 

In order to graphically illustrate the behavior of the dif- 
ferent particular models, we first introduce the dimensionless 
surface density of the disks, defined as 



2 37rMG -2 

---g^^^(4^3^^ 
^^2^ 105^^2(5^4 _ ^2^2 



M - 

'Sm{R) = Tr'^miR), 



(29) 
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for < -R < 1. So, in Fig. [T] are depicted the dimensionless 
surface mass densities for the models corresponding to 
m = 1, 8. As we can see, the disks with higher values of 
m present a mass distribution that is more concentrated at 
the center and less at the edge. Accordingly, these disks can 
then be considered as appropriated models of galaxies with 
a central bulb. 

Now, in order to graphically illustrate the behavior of 
the circular velocities or rotation curves, we introduce the 
dimensionless quantity 



Vm{R) = \J^Vm(R). (30) 

We plot, in Fig. [5] the dimensionless rotation curves for the 
models corresponding to m = 1, 10. The circular velcotiy 
corresponding to m = 1 is proportional to the radius, repre- 
senting thus a uniformly rotating disk. On the other hand, 
for m > 1, the circular velocity increases from a value of 
zero at the center of the disks, until it attain a maximum 
at a critical radius and then decreases to a finite value at 
the edge of the disk. Also we can see that the value of the 
critical radius decreases as the value of m increases. 



5 CONCLUDING REMARKS 

We presented an infinite family of axially symmetric thin 
disks of finite radius obtained by means of a particulariza- 
tion of the Hunter method. The disk models so obtained are 
generalizations of the well known Kalnajs disk, which cor- 
responds to the first member of the family. The particular- 
ization of the Hunter model was obtained by requiring that 
the surface density was a monotonically decreasing function 
of the radius, with a maximum at the center of the disk 
and vanishing at the edge, in such a way that the mass dis- 
tribution of the higher members of the family were more 
concentrated at the center. 

We also analyzed the rotation curves of the models and 
we find for the first member of the family, the Kalnajs disk, 
a circular velocity proportional to the radius, representing 
thus a uniformly rotating disk, whereas for the other mem- 
bers of the family the circular velocity increases from a value 
of zero at the center of the disks until reach a maximum at 
a critical radius and then decreases to a finite value at the 
edge of the disk. Also we find that the value of the critical 
radius decreases as the value of m increases. 

We believe that the obtained thin disk models have 
some remarkable properties and so they can be considered 
as appropriated realistic flat galaxy models, in particular if 
the superposi tion of these thin disks with appropriated halo 
distributions l|Binnev fc Tremaind l|l98i1 )) is considered. We 
are now considering some research in this direction. We are 
now also working in the non axially symmetric generaliza- 
tion of the here presented disks models and also in the obten- 
tion of the relativistic generalization of them for the axially 
symmetric case. 
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Figure 1. Surface density Em as function of R for generalized 
Kalnajs disks models with m = 1 (bottom curve) until m = 8 
(upper curve). 




0.2 0.4 0.6 0.8 1 



Figure 2. Rotation curves Vm as functions of R for generalized 
Kalnajs disks models with m = 1 (right line) until m = 10 (upper 
curve). 
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